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Particle  Trajectories  and  Potentials  in  a 
Plane  Sheath  Moving  in  a  Magnetoplasma 


1.  INTRODUCTION 

In  general,  an  analytic  self  consistent  solution  cannot  be  obtained  for  the 

interaction  of  a  charged  body  with  a  magnetoplasma.  The  only  exception  found,  the 

case  treated  here,  is  that  of  an  infinite  planar  sheath  parallel  to  the  magnetic  field. 

1  2 

(For  an  infinite  cylindrical  magnetron.  Page  and  Adams  ’  obtained  the  cutoff 

radius,  magnetic  field,  and  a  good  approximation  to  the  trajectories  and  electric 

field.  )  The  solution  for  a  crossed-field  gap  (plane  diode,  planar  magnetron)  is  old 
3  4 

(Benham,  Birdsall  and  Bridges  and  references  therein).  Instead  of  emission  from 
one  plate  toward  a  parallel  plate,  here  only  one  plate  (object  plane)  moving  relative 
to  a  plasma  is  considered.  Nevertheless,  for  the  model  considered,  the  cases  are 
close  enough  for  the  solutions  to  be  basically  the  same.  For  simplicity,  we  assume 
the  sheath  edge  to  be  sharp  and  parallel  to  the  object  plane,  collisions  to  be  negligible. 


(Received  for  publication  6  July  1987) 
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ami  all  <>l'  thi'  attracted  particles  to  have  the  same  sheath  entry  angle  and  perpen- 
djrnlar  component  of  energy.  Whether  or  not  the  sheath  edge  may  he  considered 
sharp  depends  on  the  relative  values  of  various  scale  lengths  such  as  Debye  length, 
repelled  species  ambient  gyro  radius,  and  total  sheath  thickness.  By  these  criteria, 
tor  the  ionospheric  applications  envisioned,  the  sheath  edge  is  not  sharp.  Even  as 
a  poor  approximation,  however,  the  physical  trends  of  the  results  obtained  here 
art*  exported  to  be  correct.  Without  motion  between  the  plasma  and  object,  and  the 
collision  mean  tree  path  much  larger  than  the  gvroradius  as  in  the  K  layer  of  the 
ionosphere,  the  steady  state  current  density  to  such  a  surface  would  be  negligible. 
This  is  due  to  lack  of  a  mechanism  for  charged  particles  continually  to  come 
within  two  gyro  radii  of  the  sheath  edge  where  their  gyration  could  cause  impact. 

(A  finite  plane  would  actually  receive  appreciable  current  along  the  magnetic  field 
lines  to  and  near  the  edge  of  the  plane.  )  W  ith  motion,  the  solution  is  slightly 
different  from  that  for  a  planar  magnetron.  Here,  cutoff  or  the  amount  of  current 
turned  back  is  related  to  a  function  of  parameters  called  q  as  defined  and  used  by 
Hinson,  J  for  example.  The  range  of  incidence  angles  for  a  moving  object  is  calcu¬ 
lated,  the  stability  of  a  cutoff  sheath  is  examined,  and  conclusions  are  drawn. 


5.  I.inson,  R.  H.  (1969)  Current-voltage  characteristics  of  an  electron-emitting 
satellite  in  the  ionosphere,  J.  Geophys.  Res.,  74:2368-2375. 
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2.  TRAJECTORIES  AND  FIELDS 

Figure  1  shows  the  coordinate  system 
chosen.  The  sheath  edge  is  at  yQ,  y*. 
is  perpendicular  to  y,  y*,  and  is  assumed 
to  be  sharp:  the  plasma  is  assumed 
neutral  with  no  electric  field  for  y  =  y  , 

j  jq 

and  the  repelled  species  is  assumed  not 
to  penetrate  to  y  >  y  .  E  ,  the  electric 
field  at  y  -  y  ,  is  carried  through  the 
analysis  anyway  for  generality  and  for  the 
stability  treatment,  is  the  component 
perpendicular  to  B  of  attracted  species 
thermal  velocity  in  the  plasma  frame. 

The  plasma  moves  toward  the  object  with 
a  relative  velocity  v^  along  y.  vq  =  v^+  v^. 
The  attracted  species  is  assumed  to 
enter  the  sheath  at  an  angle  a*  in  the 
plasma  frame  and  a  in  the  object  frame. 
Calculation  is  made  in  the  object  frame. 
Figures  2  and  3  show  qualitatively 
examples  of  trajectories;  near  the  sheath 
edge  for  B  =  zB  and  B  >  0,  positive 
particles  move  in  the  counterclockwise 
and  negative  in  the  clockwise  sense.  A 
dot  over  a  quantity  means  its  time 
derivative.  The  electric  field  is  in  the 
y  direction,  and  the  cylindrical  radius, 

R,  is  perpendicular  to  B.  So, 


Z,  ±  B 


Figure  1.  Coordinate  System 
(a)  Plasma  Frame  and  (b)  Object 
F  ramp 


Figure  2.  Penetrating  Trajectories 


Figure  3.  Reflecting  Trajectories 


R  =  xx  +  yy,  B  =  pzB,  p  •-  ±  1,  B  >  0,  E  =  yE  =  -yd  <p/dv  . 


The  initial,  thermal,  drift  velocities,  and  entranc  e  angle  arc-  related  as  follows: 


x  -  v  cos  o:  =  v  cos  (v*.  v  =  v  sin  a  =  v  sin  at*  +  v 
o  o  g  •  o  o  g  d 


9  «9  2  9  2 

v  --x  4-  v  ,  v  -  x  +  (y  -  v  , )  =  v  -2vv.sincv+v, 

o  O  '  o  g  o  o  d  o  o  d  d 


The  differential  equations  for  the  attracted  particles  are  the  motion,  c  n<-rgv,  and 
Poisson  equations: 


R  =  S  —  [E  +  R  X  B] 
m 1  ‘ 


•  2  2e  2  2e 

R  +  s  — —  o  -  v  +  s —  m  -  const, 
m  r  o  m  ro 


7*  “  • 

dy  o 


where  s  =  ■+  1  and  -1  for  the  attracted  species  being  ions  (assumed  positive-,  singly 
ionized)  and  electrons  respectively,  e  is  the  magnitude  of  the  electronic  charge,  and 


m  is  the  mass  of  the  attracted  species.  Except  for  <  ,  the  subscript  o  denotes  the 


mitial  value,  the  value  at  the  sheath  edge.  The  time  t  is  chosen  zero  when  the 
particle  enters  the  sheath  edge.  Bet 


eB/m,  9  =  oil,  v  r  oir,  X  =  x/r,  Y  =  v/r, 
g  g  g  g 


dX  •,  .  _  dY  •,  *  _  2eo 

x  3  X,V  '  3  JcTy/v  s’  *  3 

K 


d  i  T.  2c 
dY  m  <jl)  v 

g 


oi  “  n  m  .  r  r  -i2 

q  — J?-  =  — - — it  =  t;  t- -  ,  M  =-  V  ,/v 

oT  t  ir  2  LadJ  d 

o 


V  =  V  /v 
o  o  g 


2  2 
X  1  f  Y  1 
o  o 


1 1  +  2IVI  sin  c»*  +  1VI 


X  1  =  cos  a*  -  V  cos  cv,  Y'  -  sin  a*  4  M  -  \  sin  cv 
o  o  o  o 


Differentiation  by  0  is  denoted  by  a  prime;  denotes  the  ambient  plasma  number 
density;  VI  is  the  Mach  number;  A  is  the  Debye  length.  Normalizing  thus,  the 
differential  equations  become 

\"  =  SpV  ' 

(2M 


v "  c  Ik  -  sp.v 


2  2^  2  . 
V  4  '  +  4>  -  1  +  2M  sin  a/*  +  M  + 


(3) 


■> 

d'i. 


dV  ‘ 


-  2q  — 
M  n 


(4) 


To  treat  both  species  simultaneously,  we  shall  take  p  =  s.  Then  Eqs.  ( 1 ' )  and  (2M 
become 


X"  =  V 


(1) 


V"  = 


X' 


(2) 


q  >  0,  4>  -  $  5  0,  and  E  i  0.  Thus  attracted  electrons  have  the  same  normalized 

o 


trajectory  as  attracted  ions  if  B  has  (and  <p  and  E  have)  the  opposite  sign. 
Integrating  Eq.  (1), 


X' 


cos  (»*  +  Y  -  Y 


(5) 


Continuity  of  current  provides  the  additional  expression  needed  to  determine 
the  density.  This  expression,  two  of  the  first  three  equations,  and  Eq.  (4)  then 
constitute  a  close  set.  First,  we  consider  three  possible  cases  concerning 
tra  jectories.  (See  Section  3.  ) 


(A) 


(H) 

(C) 


q  is  sufficiently  large;  Y(0  )  increases  and  never  decreases  with  0  ,  as 
in  Figure  2.  All  of  the  current  is  transmitted  and  is  assumed  to  be 
collected. 

Between  the  two  limiting  values  of  q,  a  fraction  of  the  particles  return 
and  the  rest  are  transmitted. 

q  is  sufficiently  small,  and  object  |$|  and  thus  sheath  edge  to  object 
distance,  are  sufficiently  large;  the  particles  all  return  as  in  Figure  3. 
This  is  the  cutoff  case  with  no  net  current. 
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How  do  the  return  and  forward  velocities  compare?  Since  -  4-  i 'i  I  and,  from 
Eq.  (5),  X'  =  X'(Y),  then  from  Eq.  EH,  for  the  same  value  of  N  ,  y  ' .  -V  ,  where 
>  0  and  Y'  <  0.  So,  the  return  part  of  a  trajectory  is  symmetric  to  Hie 
Y'  >0  part. 

Since  the  y  component  of  incident  particle  flux  equals  the  sweepout  rate  ami 
the  steady  state  current  is  conserved. 


n  Y '  v  =  n  Y '  v  -  n  v  .  , 

+  +g  O+CH-g  d 


where  n+  and  n_  denote  the  number  density  of  particles  having  Y'  '>  0  and  Y'  <  0 
respectively.  This  assumes  that  any  returning  particles  exiting  the  sheath  edge  do 
not  reenter  the  sheath.  Then  n  =  Mn  /Y1  .  This  Y'  >  0  part  of  the  current 

0+  “O  o 

density  is  constant  up  to  a  turning  point,  0  (V  =  0)  a  r  ,  if  a  turning  point  exists. 
Let  F  -  1  be  the  fraction  of  the  incident  current  density  that  is  returned  if  the 
collector  is  beyond  a  turning  point.  Then  the  return  particle  current  density  at  the 
same  value  of  Y  is 


-n  Y'v  -  n  Y'v  =  (F  -  l)n  v, 
*  -  g  -  +  g  d 


Since  steady  state  is  assumed  and  no  charge  sources  or  sinks  exist  in  the  sheath, 
the  net  y  component  of  flux  is  the  same  throughout  the  sheath.  This  then  is 
[1  -  (F  -  l)ln  v  ,.  F  =  l  for  Case  A,  1  <  F  <  2  for  case  B,  and  2  for  case  C. 

Eq.  (7)  times  n+  divided  by  Eq.  (6)  shows  that,  for  the  same  value  of  Y, 

n  =  (F  -  l)n  .  (8 

+ 

Before  reaching  any  turning  point, 

n  =  n  +  n,  -  FMn  /Y'  ,  (9 

from  the  sum  of  Eqs.  (6)  and  (7)  divided  by  Y_J_v  .  So,  for  all  three  cases,  for 
Y'  >  0,  Eq.  (4)  becomes 


d  $  2FMq 

- 2 - V1 — 

JU  **  1 


Multiplying  by  dY  =  Y'd0  and  integrating; 


6 


-  d$>/dY  =  E  =  Eq  +  2FMq9 


(11) 


From  Eqs.  (2),  (5),  and  (11), 

Y"  =  FMq0  +  a  +  Y  -  Y  ,  (12) 

o 

where 

a  =  E  / 2  -  cos  a*  .  (12) 

o 

This  quantity  has  two  parts:  that  due  to  the  non-space-charge-limited  initial 
electric  field,  and  that  due  to  the  angle  of  the  particle's  incident  velocity  with  the 
sheath  normal. 

The  Ansatz  Y  =  -A  cos  (0  +  0)  +  ('.0  +  1J  is  substituted  into  Eq.  (12).  This 
and  the  initial  conditions  lead  to  the  solution  rO"  Y'  ?  0. 

Y  ^  Y  f  a  -  A  cos(0  +  8  )  +  FMq9  (14) 

o 


Y  1  =  A  sin  (9  +  )  +  FMq 


( i  r. ) 


Y"  =  A  cos  (0  +  (i  ) 


(Hi) 


A  = 


■ya2  +  (YJ  -  FMq)2 


cos  -  a/  A  - 


4- 


( 


Sin  a  * +  M  -  F  M  q 

E  /  2  -  cos  a* 
o 


-1/2 


(IV) 


(  1C) 


sin  -  !•  M q)/A  .  (I'D 

The  quantity  a  is  seen  to  be  the  Y  component  of  the  initial  normalized  u>  <  deration. 

4 

Eq.  (14)  agrees  with  Birdsall  and  Bridges,  |Section  5.04.  Eq.  CO],  up  to  the 
coefficients  of  the  functions  of  0.  Substituting  Eq.  (11)  into  Eq.  (a),  integrating, 
and  using  the  initial  conditions  leads  to 

X  =  X  +  V  '  -  FMq  -  A  siri(0  <  d)  +  <E  1 2)0  +  (!•  \lq/2)02  •  20) 

o  o  o 


X'  -  -  \  cos(0  1  d  )  +  E  /2  +  I'  M qu 

o 


i  2  1 ) 


X 


\  sin(9  1  /(  )  *  1  \;<)  . 


<  2  2) 


cuuv«r w 


i  n*  t*»,i|i,i 


From  Eq.  (3) 


r2  _  vi2  _  V  i2 


<J>  =  $  +  1  +  2M  sin  <y*  +  M  -  X'  -  Y ' 
o 


With  M  =  0,  this  agrees  with  Birdsall  and  Bridges  [Section  5.05  Eq.  (7)]. 

Forty*  <  7)  12,  X'  >  0  for  all  values  of  Y.  Forty*  >77/2,  X'  <  0  for 
Y  -  Y  <  -cos  a*  and  X'  >  0  for  Y  -  Y  >  -cos  a*.  For  O'*  =  77/2,  the  trajectory 


is  nearly  straight  for  large  q,  wave  shaped  for  q^  <  q  ~  q  (see  Sections  3  and  4), 


U-shaped  for  reflected  particles,  and  less  than  half  a  U  if  collected  before  reflection. 

In  summary,  the  trajectory  of  each  particle  is  given  by  Eqs.  (20)  and  (14)  for 
Y'  =  0  and  is  symmetric  for  V  <0  (the  particles  turned  back);  the  electric  potential 
and  field,  before  a  turning  point  if  any,  are  given  by  Eqs.  (23)  [with  Eq.  (21)  and 
Eq.  (15)]  and  Eq.  (11)  respectively;  for  case  C,  the  electric  field  beyond  the  turning 
point  equals  the  value  at  the  turning  point. 


3.  THE  TURNING  POINT 


If  the  attracted  particles  can  come  to  a  point  where  Y'  =  0,  called  "the  turning 
point",  generally  the  magnetic  field  reflects  at  least  some  of  them  back  to  the  sheath 
edge.  The  limit  of  q  for  this  occurrence  and  the  value  of  6  ,  X,  Y,  and  at  the 
turning  point' a  re  now  calculated. 

From  Eq.  (15),  case  A  (F  =  1)  occurs  if  and  only  if  A  s  Mq.  From  this  and 
Eq.  (17),  the  minimum  value  of  q  for  case  A  is 


22—?  —  2 
a  +  Y'  E  /4  -  E  cos  a*  +  V 
o  o  o  o 


(sin  q*  +  M) 


q^  is  a  function  only  of  initial  (sheath  edge)  electric  field  (if  any),  entrance  angle. 


and  M. 


For  q  1  q  .  from  Eq.  (15),  T  ,  the  value  of  9  which  makes  Y'  first  vanish, 
m 


satisfies 


sin  (T  +  (3  ) 


-F  Mq 


COS  (T  +  /3  )  = 


(FMq)' 


® 


where 


/a2  -  (FMq)2 


Y  '(Y  '  -  2F  Mq' 
o  o  ^ 


The  sign  of  the  right-hand  side  of  Eq.  (26)  is  due  to  T  +  fS  h  iving  to  be  in  the  third 
quadrant.  Expanding  sin  r  =  sin[(T  +  (3  )  -  /3  ],  cost  =  eos|(T  +  (i )  -  (3],  and 
using  Eqs.  (18),  (19),  (25),  and  (26)  leads  to 


sin  T  = 


S(Y^  -  FMq)  -  aFMq 


cos  T  =  - 


aS  +  FMq(Y^  -  FMq) 


By  setting  9  =  t  and  substituting  Eqs.  (25)  and  (26)  into  Eqs.  (20)  and  (14),  we  find 


the  coordinates  of  the  first  stationary  point  for  q  §  q^  to  be 


X  =  X  +  Y  '  +  T  +  £^2  r  2  , 
r  o  o  2  2 


Y  -  Y  +  a  +  S  +  FMqr  . 
r  o 


This  is  usually  the  reflection  point  of  the  particle.  The  above  results  lead  to  the 
following  deductions  as  q  goes  from  a  large  to  a  small  value. 


Case  A: 


q>q  :  S*  <  0,  F  =  1,  A  <  Mq,  the  only  case  for  which  no  stationary  points 


exist. 


q  =  q  :  the  limit  of  case  A,  S  =  0,  F  =  1,  A  =  Mq,  sin  T  =  -a/ (Mq),  and 
n  m  _ 


V2  2 

1  -  a  /(Mq)  .  Stationary  points  also  occur  at 


9  =  T  +  2  tri,  i  =  1,  2,  3,  . . .  until  the  particles  are  collected. 


Case  B: 


qm/2  <  q  <  qm:  =  0,  F  =  qm/q,  A  “  FMq,  and  T  ,  Xr>  Yr,  A,  and  the 

particle  trajectories  are  independent  of  q.  If  the  stationary  point  is  before  the 
object,  it  is  a  reflection  point  for  F-l  fraction  of  the  particles;  the  remaining  2-F 
fraction  of  the  particles  penetrate  to  the  object. 


-VT  v/.  .s  A.N.V.  .v-.v. 
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Case  C: 


q  =  q^/2:  the  limit  of  Case  C,  the  maximum  value  of  q  for  complete 


cutoff;  all  particles  are  reflected  by  the  magnetic  field  unless  the  turning  point  is 

2 

beyond  the  object  surface.  S  =  0,  F  =  2,  A  =  2Mq. 

2 

q  <  q  /2:  S  >0,  F  =  2,  A  >  2Mq.  and  T,  X  .  and  Y  depend  on  q. 
m  r  r 

0,  what  is 


What  are  the  limits  of  o',  T  ,  q  ,  and  the  reflection  points  as  Y 

m  o 


the  range  of  a*,  and  what  is  the  corresponding  range  of  a  ? 
For  M  <  1,  the  attracted  species  will  enter  the  sheath. 


0  <  a  <  if,  for 


—  77/2  <  -sin  1  M  i  a*  5  V  +  sin  1  M  <  3;r/2. 


If  also 


-2  y 1  -  M2  <  Eo  <  2  V 1  -  M2,  then  for  0  s  sin"1  M  <  77/2, 


-1 


( 1 )  as  a*  —  -sin  M: 


a<0,  Y'-*a-‘T-*-0,  q  10 ,  X  -X  and  Y  —  Y  as  they  must,  and 
o  Mm  r  o  r  o  J 


1 


(2)  as  a*  -*  7i  +  sin  M: 


a  >  0,  Y'  -  0,  a  -»7,  )i  i  Tl  2t,  q  —  *> ,  X  *  X  +qMr  , 
o  m  r  o 


and  Y  -*  Y 
r 


o  ” 


1  -  M  +  2qMT 


Physically,  the  magnetic  field  swings  the  particles  away  from  the  sheath  in 
(1)  and  into  the  sheath  in  (2). 

For  M  >  1: 


0  <  cos”  1  1  /M  s  a  =  V  -  cos  1  1  /M  <  77  , 


1 


when 


—77  / 2  +  cos”1  1/M  g  a*  g  -tj/2  -  cos”1  1/M  and 


and  when 


-77/2  +  cos  1  1/M  5  a*  ^  3jt/2  -  cos”1  1/M  . 


10 


ZKA 

X'K 

< 


For  E^  -  0  [space  charge  limited) 


V 


1  -  M 


m  2M  sin  ft  2M  (sin  a*  +  VI) 


+  1 


(32) 


This,  for  the  three  limits  of  M,  is  as  follows. 


( 1 )  For  M  «  1,  q.„  =  1  /(2MV '  )  +  1.  Fora  given  value  of  M,  the  minim 

1  n  1  ° 

value  is  q  =  tj  -rs  +  1  for  «*=(>  =  7) 72.  The  maximum  value  is  *>for 
m  2  j 

rv*  =  -sin  *  M,  <>  =  0,  or  tt. 


um 


(2)  For  M  =  i  (strictly,  for  |  M  -  1 1  «  1  +  sin  a*):  q  -  1.  For 

m 

sin  a*  4  -1  and  M  =  1,  q  =  1. 

m 

(3)  For  M  »  1,  qm  =  1/2. 

F  rom  Eq.  (32),  three  qualitatively  different  types  of  curves  of  q  (M)  occur: 

^  .7  r-  ^m 

Type  1:  sin  re*  =  -1,  a  singular  case,  v  =  -v  y,  shown  as  curve  1  of  Figure  4. 

g  • 

The  thermal  velocity  is  equal  but  opposite  to  the  drift  velocity,  so  y  remains  zero. 


Tvpe  2:  -1  <  sin  a*  <  0,  v  <0,  between  the  limit  of  curve  2  as 

gy 

sin  a*  —  -1  (the  M  =  1  line  and  curve  1)  and  curve  3  of  Figure  4.  In  the  plasma 

frame,  the  particle  is  moving  away  from  the  sheath;  however,  the  sheath  is  moving 

faster,  so  it  catches  up  to  the  particle.  A  minimum  occurs  at  M  =  Mm  >  1  and 

q  =  q  ,  where 
m  mm 


sin  a'-  =  - 


M  +  1 
m 


id  --Vi 


A  maximum  occurs  at  M  =  <  1,  which  also  satisfies  these  two  equations.  But 

q  (Mm)  <  0;  which  is  physically  meaningless  except  that  the  sheath  can  never 
catch  up  to  the  particle. 

Type  3:  sin  a-'  =  0,  v  5  0,  between  curves  3  and  4  of  Figure  4. 

gy 

What  does  this  tell  us  about  the  effect  of  varying  v^  from  0  (M  -  °°)  to  the 
maximum  value  (minimum  M  -  the  larger  of  0  and  -sin  o*)? 

(1)  For  type  1 

(a)  1  <  q;  case  A.  (all  incident  particles  collected) 

(bl  0.  5  <  q  <  1:  case  B  (some  particles  collected)  goes  to  case  A. 

(c)  q  <  0.  5:  case  C  goes  to  case  A.  Case  A  exists  up  to  M  =  1. 

If  M  <  1,  then  q^  <  0,  which  is  physically  meaningless  except 
that  the  particles  can  never  cross  the  sheath  edge. 

(2)  For  type  2 


(a)  0.  5  <  q:  case  A  goes  to  case  B  at  M  <  M 


q  =  q,  then 
m 


finally  to  case  C  (no  particles  collected). 

(b)  0. 25  <  q  <  0.5:  case  B. 

(c)  q  <  0.  25:  case  C  may  go  to  case  B. 

Either  of  the  latter  two  may  go  to  case  A  at  M  <  M  and  back  to  case  B  at  M  <  M 

m  m 

before  finally  going  to  case  C.  For  all  of  type  2,  this  final  case  C  separates 

farther  from  case  B  without  limit.  That  is,  if  E  =  0, 

- 1  ° 

0  <  sin  M  <77/2,  and  -sin  O'*  <  M  -*■  -sin  dr*,  then  (i)  q  -*«*>,  (ii)  X  "X  and 

- 1  2  m  r  0 

Y  -*■  Y  for-sin  —a*,  and  (iii)  X  —  X  +  qMr  and 
r  o  _  r  o 

V2  -1 

1  -  M  +  2qMr  ,7/  <  r  <2.7  for  v  +  sin  M  —  O'*.  So,  for  sufficiently 
large  yg  -  yQ  (thus  magnitude  of  object  potential),  finally  the  particles  are  reflected 


at  y  =  y. 


If  M  <  -sin  (V*,  then  q  <0,  which  is  physically  meaningless  except  that 

m 


the  particles  move  away  from  the  sheath  ■edge,  and  so  can  never  enter  the  sheath. 


« *4 ' i'i  *>.*S  !*>  rtWlVf 


avXuonTXinn 


(3)  For  type  3 

(a)  0.  5  <  q:  case  A  goes  to  case  B,  then  finally  to  case 

(b)  0.  25  <  q  <  0.  5:  case  B  goes  finally  to  case  C. 

(c)  q  <  0.  25:  case  C  remains  Case  C.  Case  C  separates  farther  from 

case  B  without  limit:  q  as  M  -*•0.  At  first,  this  may 

m 

contradict  intuition  since  the  effect  of  magnetic  field  must  decrease. 

We  see  from  Eq.  (31),  however,  that  y^  approaches  a  function 

that  increases  linearly  with  as  M  -  0;  thus  the  particles 

penetrate  the  sheath  without  limit. 

If  F  =  1,  E  =0,  and  q  =  q  ,  then 
o  m 

T  =  2  a  .  (35 


4.  RANGE  OF  INCIDENT  ANGLES 

To  get  the  range  of  incident  angles, 
we  consider  the  gyration  phase  angle 
before  the  particle  enters  the  sheath  at 
time  t  =  0,  taken  to  be  when  the  circle  of 
gyration  becomes  tangent  to  the  sheath 
edge.  As  before,  a  sharp  sheath  edge 
is  assumed,  and  the  plane  surface  and 
plasma  are  coming  together  along  y 
with  a  relative  velocity  v^.  The 
coordinate  system  is  chosen  in  the 
object  frame  with  initial  conditions  as 
in  Figure  5:  the  origin  on  the  sheath 
edge  and  the  gyrocenter  a  gyro  radius, 
r  .  before  it.  The  perpendicular 
component  of  velocity  in  the  plasma 
frame  is  v  =  wr  .  In  general,  this  /3 
is  different  from  the  p  used  previously. 
The  particle  then  is  at 


\  .  v« 

\  / 9 


sheath  V 


/ 

edge 


object  surface 


x  =  r  sin  (cut  +  /3  ) 


Figure  5.  Coordinate  System  and 
Initial  Conditions 


y  =  v^t  -  rg(l  "  cos  (cut  +  0 )]  . 
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(37') 


(36) 


Normalizing  as  in  Section  2, 


X  =  sin  (8  +  ) 


Y  =  Mfl  -  1  +  cos  (8  +  P  )  . 


Let  the  subscript  o  designate  the  value  for  which  Y  =  0,  and  a  prime  designate  the 
derivative  with  respect  to  8  .  The  entrance  angle  is  determined  from  the  values 
at  the  sheath  edge  of  the  derivatives  of  Eqs.  (37)  and  (36): 


tan  g 


M  -  sin  (8  f  /3  ) 
cos  (8 o  +  (3 ) 


Case  1.  M  §  1. 


Let  /3  be  the  value  of  (i  which  makes  8  Q  the  smallest  non-negative  value,  &  Qy 

for  which  Y  1  =  0,  and  let  8  „  be  the  next  larger  value  of  8  with  B  =  B  . .  The 
o  o2  a  o  1 

extrema  of  o  occur  for  (3  =  (3  :  the  minimum,  0,  for  8  =  0^,  anc*  the  maximum 
for  8  --  8  9.  Equating  the  numerator  of  Eq.  (38)  to  zero. 


sin  (8  4/3^=  M  . 


Using  this  in  Eq.  (37)  =  0, 


8  Ql  =  (1/M)  [  1  -  cos(fl  +  (3  )]  =  M / [  1  +  Y 1  -  M  ] 


M 

% 


From^1  =  (flQl  +  0j>  -eQl. 

(3  =  sin"1  M  -  M/[  1  +  -^1  -  M2] 


A  particle  with  a  slightly  larger  value  of  (3  must  gyrate  through  an  additional  angle 
6  o 2  ~  ®  0\  dess  than  a  gyration)  before  reaching  the  sheath  edge,  and  then  enters 
it  at  q  jyj,  the  largest  value  of  a. 

8  q2  =  (1/M)[1  -  cos(9o2  +  0  j)]  .  (42) 

So,  the  range  of  a  is 


wmmssm 


j  M  -  sin  (0  o2  +  /3  ^ 


cos  (0  o2  +  V  jT 


As  M  -  0:  0  ol  */3  1  -*  M/2,  0  o2  -*  2  n  -  f>M,  and  nr  -  2\'j;  M.  At 


M  =  1:  ft  ,  -  9  „=1,  /3=7T/2-l,  and  0!  =  i . 

ol  o2  1  M 


Case  2.  M  =  1 . 


Only  one  value  of  9  exists  for  each  value  of  .  Solving  10 q.  (37) 


0  for 


cos  (9  +  /3  )  and  substituting  into  Eq.  (38), 


fM9  (2  -  M0  ) 


O  =  tan 


1  -  Mf 


,0^0  s  2/M 
o 


Differentiating  this  with  respect  to  9  o>  equating  to  0,  solving  for  0  ,  and  sub¬ 


stituting  into  Eq.  (44)  leads  to  the  minimum  and  maximum  values: 


I M  -  1  i  (O  I  -  tan 


/ M"  -  1  . 


As  a  numerical  example,  consider  a  satellite  in  a  low  earth  orbit  and  assume 


that  T  =  1000°  K,  v  .  =  7  km/sec,  n  =  101 1  m  ^  and  B  =  0.  5  X  10  ^  T.  Then,  for 

d  30 


electrons,  M  =  0.05,  0  A  rt  5  45°,  14  <  q  <  «o ,  and  q  -  4,1.  For  O  ions, 
c  rn  c 

MQ+  =  8.4,  83°  =  a  =  97°,  0.44  <  q^  <0.56,  and  q  =  1.  2  X  10  .  So,  the  theory 


predicts  that  electrons  with  the  mean  thermal  velocity  would  be  reflected  unless 
the  distance  from  the  sheath  edge  to  the  object  surface  were  less  than  the  reflection 
distance,  that  is,  unless  the  magnitude  of  the  object  potential  were  small,  where¬ 
as  ions  would  be  practically  unaffected  by  the  magnetic  field.  For  attracted  elec¬ 
trons  with  the  numerical  values  above.  Figure  6  shows  normalized  trajectories 
from  the  origin  and  Figure  7  normalized  sheath  potentials  for  three  incident  angles, 
all  case  C.  The  "trajectory"  for  a  =  0,  curve  0  in  Table  1,  is  the  point  at  the 


origin. 


m 

is 


B 


II 


VVv 


K* 


Table  1.  Parameters  of  Figures  6  and  7 
M  =  0.  05,  q  =  4. 1 


<y* 

a 

% 

P 

T 

X 

r 

Y 

r 

C  u  rve 

0 

-2.  9° 

0° 

QO 

202.3 

0° 

0 

0 

1 

7.  2 

10 

53.  0 

193.3 

10.4 

0.  18 

0.  016 

2 

27.  5 

30 

20.5 

173.4 

33.  9 

0.  58 

0.  15 

3 

43.  0 

45 

14.  6 

156.  3 

54.  6 

0.  92 

0.  35 

The  x's  in  Figure  7  are  at  Y  =  V  .  For  Y  >  Y  ,  the  curves  are  linear  since  there 
is  no  charge  and  E  is  constant.  Figures  8  and  9  are  trajectories  and  potentials 
respectively  for  0+  ions,  case  A.  The  curves  for  the  largest  value  of  q  are  with 
the  parameters  in  the  first  part  of  this  paragraph,  and  are  replotted  in  Figures  10, 
11,  and  12.  With  no  magnetic  field.  Figure  10  and  curve  2  of  Figure  11  would  be 
straight  vertical  lines.  As  expected,  they  show  that  the  geomagnetic  field  deflects 
the  ions  very  little.  Also  as  expected.  Figures  11  and  12  show  a  narrow  range 
of  incidence  angles  and  potential  distributions  for  0+  ions.  Although  this  calculation 
is  not  realistic  since  the  sheath  edge  is  not  really  sharp  and  there  is  a  distribution 
of  velocities,  it  shows  that,  for  low  earth  orbit  conditions,  the  predominant  values 
of  a  are  low  for  electrons  (object  positively  charged)  and  near  90°  for  ions  (object 
negatively  charged).  Also,  as  expected,  generally  ions  are  affected  very  little 
whereas  electrons  are  greatly  affected  by  the  magnetic  field. 
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5.  STABILITY  OF  A  CUTOFF  SHEATH 


Section  3  shows  that,  for  this  simpli¬ 


fied  case  and  q  <  q  /2,  all  of  the 
m 


attracted  particles  are  reflected  away 
from  the  object  plane  at  a  certain 


Z.B 


y  =  y„  (if  y  >y  ,  see  Figure  13), 

I  o  r 


making  a  vacuum  insulation  gap.  The 
theory  also  shows  that  the  difference  in 
object  and  reflection  point  potentials. 


E  =  0 


<ps  -  is  proportional  to  y  -  y^; 


thus  the  gap  increases  with  increasing 


sheath  edge 

1 - Tr. 


This  is  consistent  with  the 


repelled  species  being  repelled  farther 


CHARGE 


from  the  object  for  <p  larger. 


Is  this  result  of  the  theory  correct? 
Here,  we  consider  this  question,  h  irst, 
we  assume  a  gap  and  determine  whether 
or  not  it  tends  to  close  (decrease  in 


gap  yr 

object  surface 


thickness).  Let  be  the  surface  charge 


density  on  the  object  and  be  the 


sheath  charge  in  a  column  of  unit  cross 


section  area  from  y  to  v. 

■'  o 


Figure  13.  Insulation  Gap 


8  FMn 


a  =  /  snedy  =  se  / 

c  y  o 

'  o 


—  v  dt 


sFn  ev,  9  loo, 

jo  q 


from  Eq.  (9)  and  the  normalization  in  Section  2.  For  y  %  v  <  v  : 

r  S 


a  =  cr  .  =  sFn  ev,  t/u>  , 

c  CM  30  d 


the  total  charge  per  unit  area  in  the  sheath.  So,  in  general,  a  -  cr^  and 


$  "*  $  as  Y  “*■  Y  .  We  start  with  cr  -  -cr  „„  (10  =  0)  and  the  object  isolated, 

a  o  o  s  s  c.M  o 

A  variation  of  Y  ,  6Y  ,  with  Y  fixed  causes  no  variation  in  cr  while  v  '  \  and 
o  o  s  e  M  '  s  '  r 


no  variation  in  since  the  current  to  the  object  is  zero  (although  the  object 


potential  varies).  (The  same  results  would  be  obtained  by  varying  y  instead  of  y 


Therefore,  5  Eq  =  0  so  that  6a  =  5  A  =  6S  =  0.  Also,  from  Eq  (47),  6  r  -  0.  The 


maximum  magnitude  of  magnetic  field  due  to  the  current  is  at  \  -  y  ^  or  v  and  may 
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be  shown  to  be  p.  1/2,  I  =  /  j  dy,  j  -  nex.  From  Eqs.  (9),  (21),  and  the  normaliza- 

V 


U 

tion,  then  applying  the  result  to  the  case  in  Table  1  giving  the  largest  field,  that  is, 
curve  3,  we  obtain  2  X  10  ^B.  Thus,  the  magnetic  energy  is  quite  accurately  only 
that  due  to  the  ambient  magnetic  field,  so,  effectively,  it  does  not  vary  with  yQ 
or  anything  else.  The  electrical  energy  is  that  in  the  charge,  Wc»  plus  that  in  the 
insulation  gap,  W  .  Letting  a  line  over  a  quantity  designate  its  normalized  value: 


Y  Y 

— ~  I'  s 

W  =  W  +  W  ,  2W  =  f  E2dY,  2W  =  f  E2dY 
c  S  c  yJ  c  g  Y  6 

o  r 


!•'  -  2,  so. 


E  =  E  +  4  Mq  d  ,  Y  <  Y  S  Y 
c  o  o  -  ~  r 


E  =  s - ^ —  ,  Y  S  Y 

g  t  mu>v  r 

«  o  g 


-  t-2 
2W  =  f  E  Y<  d  6  . 

p  »  r* 


c  '  c 


From  Eqs.  (15),  (44),  and  (46),  W  is  seen  not  to  vary  with  Y  .  From  Eq.  (31), 

c  o  _ 

we  note  that  the  sheath  charge  thickness,  Y  -  Y  ,  does  not  vary  with  Y  .  E 

r  o  o  g 

does  not  vary  with  Y  ,  so 


6  2W  =  E20(Y  -  Y  )  =  -  E25  Y 

g  s  r  go 


from  Eq.  (31)  since  Y^  is  fixed.  Since  charges  tend  to  move  to  minimize  the 


electrical  energy,  in  accord  with  the  plates  of  a  condenser  tending  to  collapse, 

5'  tends  to  increase.  Nevertheless,  the  theory  here  indicates  that,  because  of  the 
o 

magnetic  as  well  as  the  electric  field,  the  particles  gyrate  and  drift  along  the 
x  axis  so  that,  for  q  <  q^/2and  |<pj  large  enough  to  make  Y ^  <  Y^;  an  insulation 
gap,  Y  -  Y  ,  is  maintained.  However,  the  theory  assumes  that  the  reflected 
attracted  and  repelled  particles  disappear  after  exiting  or  bouncing  off  the  sheath 
edge.  Assuming  a  plane,  sharp  sheath  edge,  calculating  two  ion  reflections  and 
three  electron  gyrations  for  the  case  of  curve  3,  Section  4  leads  to  the  following  in¬ 
ductions.  Both  species  gyrate  (in  opposite  senses),  reenter  the  sheath,  reexit,  reenter 
and  so  on.  After  first  reaching  the  sheath  edge,  due  to  the  motion  between  plasma 


and  object,  they  gain  energy  and  penetrate  more  each  gyration,  approaching  re¬ 
maining  in  the  sheath,  that  is,  with  their  minimum  value  of  y  being  y  =  0.  This 
makes  the  effective  sheath  edge  closer  to  the  object  by  the  penetration  depth  of 
the  repelled  species.  This  contradicts  the  theory  here.  Simultaneously,  both 
species  drift  (in  the  +x  direction  for  B  =  -zB).  More  importantly,  the  number 
density  of  each  species  near  the  sheath  edge  increases  by  the  ambient  number 
density  each  gyration.  So  the  electron  number  density  increases  much  faster  than 
that  of  the  ion.  If  the  object  is  positively  charged,  the  effective  q  increases  just 
as  rapidly  since  q  n  according  to  the  theory.  For  example,  for  curve  3  in 
Section  4,  the  effective  electron  q  would  be  greater  than  after  four  gyrations, 
about  4  /is.  Furthermore,  ot  increases  each  gyration,  so,  if  a  <83°,  as  it  is  for 
the  first  four  gyrations  of  the  numerical  electron  cases  of  Section  4,  q^  decreases. 
So,  case  C  "oes  rapidly  to  case  A  or  B  with  no  insulation  gap.  The  conduction 
current  will  drain  off  the  excess  electrons  so  that,  either  case  C  (or  B)  is  produced 
and  conduction  and  (partial)  insulation  alternately  occur,  or  case  B  is  maintained. 
Taking  a  distribution  of  particle  velocities  into  account,  some  of  the  higher  energy 
particles  with  sin  a*  >  0  (type  3,  see  Section  3)  will  penetrate  to  the  object.  These 
considerations  and  perhaps  others  thus  show  that  a  vacuum  insulation  gap  is  not 
expected  to  exist,  although  a  region  of  low  density  may.  We  speculate  that  essen¬ 
tially  these  phenomena  would  occur  over  most  of  a  positively  charged  plate  with 
sufficiently  large  size  and  potential. 

If  the  object  potential  and  y  are  held  constant  while  v  is  increased  (the  gap 

s  o 

decreased),  the  amount  of  charge  in  the  sheath  would  have  to  increase  (unless  K 
changes)  since  the  capacitance  per  unit  area  increases.  But  Eqs.  (28),  (29),  and 
(47)  are  independent  of  v  ,  so  this  contradicts  the  theory. 

6.  CONCLUSIONS 

An  infinitely  large,  uniformly  charged  plane  moving  perpendicular  to  its  sur¬ 
face  in  a  magnetoplasma  is  considered.  Analytic,  self  consistent  equations  of 
electric  field  in  the  sheath  and  trajectories  of  the  attracted  species  through  the 
sheath  are  obtained.  A  uniform  magnetic  field  parallel  to  the  plane  surface, 
negligible  collisions,  one  velocity  of  only  the  attracted  species  moving  into  the 
sheath,  and  a  sharp  sheath  edge  are  assumed.  Uniformity  of  magnetic  field  and 
negligible  collisions  are  generally  accurate  assumptions  for  an  object  in  low  earth 
orbit  or  higher.  Collisions  are  expected  to  have  an  appreciable  effect  if  the  object 
size  is  comparable  to  a  collision  mean  free  path.  The  assumptions  of  a  single 
velocity  and  a  sharp  sheath  edge  are  inaccurate.  Nevertheless,  the  physical  trends 
of  most  of  the  results  obtained  here  should  be  correct,  liven  though  the  object 
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moves  relative  to  the  plasma,  the  solution  is  basically  the  same  as  that  obtained 
much  earlier  for  a  planar  magnetron. 

The  effect  of  a  magnetic  field  is  determined  by  the  value  of  q  which  is  a  con- 
2 

stant  times  n^/B  .  If  it  is  much  larger  than  the  cutoff  q,  then  the  magnetic  field 
has  little  effect.  The  cutoff  q,  q  ,  is  a  function  only  of  sheath  entrance  angle 
and  M.  For  M  <  <  1,  it  goes  from  near  l/(2M)  at  vertical  incidence  to  00  at  grazing 
incidence  in  the  object  frame.  It  equals  1  for  M  =  1  and  1/2  for  M  »  1.  Thus 
q  ,  and  therefore  the  effect  of  a  magnetic  field  can  be  quite  large  at  all  values 
of  q  for  sufficiently  small  values  of  M.  This  is  expected,  since  the  effect  of  a 
magnetic  field  increases  with  decreasing  particle  energy.  For  q  <  q  and  the 
magnitude  of  the  object  potential  small  enough  for  the  reflection  distance, 
yr  -  y  ,  to  be  large  compared  to  the  sheath  thickness,  y  _  -  y  ,  the  magnetic  field 
has  little  effect. 

Numerical  examples  are  calculated  for  typical  low  earth  orbit  conditions.  For 
attracted  electrons,  the  range  of  incidence  angles  is  found  to  be  below  45°, 

M  —  0.  05,  and  q  ~  4.  The  latter  is  less  than  q  which  is  greater  than  14.  For 

.  ™  r, 

attracted  0  ions,  the  incident  angle  is  within  7  0  of  the  vertical,  M  ~8,  q  ~  10'  , 
much  larger  than  q^  which  is  0.  5.  So  the  magnetic  field  affects  electrons  greatly, 
but  it  affects  ions  negligibly. 

The  solution  presented  here  for  this  simple  case,  including  the  reflection  dis¬ 
tance,  depends  only  on  M,  sheath  entrance  angle,  and  q.  It  shows  that,  as  the 
magnitude  of  the  object  potential  increases,  the  trajectories  and  potential  function 
do  not  change;  the  sheath  just  thickens,  or,  if  there  is  an  insulation  gap,  the  gap 
increases.  However,  other  effects  (next  paragraph)  change  this  picture. 

A  stability  analysis  led  to  the  conclusion  that  an  insulation  gap  between  the 
charged  sheath  and  the  object  would  tend  to  narrow  and  thus  to  vanish.  With  a 
sufficient  magnetic  field  and  magnitude  of  potential,  however,  the  theory  indicates 
that  the  magnetic  field  swings  the  attracted  particles  clear  of  the  object  so  that  a 
gap  would  exist.  But  the  theory  assumes  that  after  the  repelled  species  hit  the 
sheath  edge  and  after  the  attracted  species  exit  the  sheath,  the  particles  disappear. 
Actually,  the  magnetic  field  causes  any  exiting  particles  to  reenter.  Assumption 
of  a  plane,  sharp  sheath  edge  and  calculation  of  a  few  gyrations  through  the  sheath 
edge  results  in  the  particles  increasing  in  number,  energy,  and  penetration  each 
gyration,  and  to  the  speculation  that  the  state  either  oscillates  between  cases  (B 
or  C)  and  A  or  remains  case  B.  Also,  the  high  energy  tail  of  a  realistic  velocity 
distribution  would  penetrate  to  the  plane.  So,  a  perfect  insulation  gap  could  never 
exist,  although  a  low  density  region  may  exist  part  of  the  time  near  the  object. 

We  speculate  that  basically  these  phenomena  would  occur  over  most  of  a  positively 
charged,  finite  plane  as  c  msidered  here,  if  its  size  and  potential  are  sufficiently 
large.  Penetration  of  the  sheath  edge  by  the  repelled  species  results  in  a 
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contradiction  to  the  theory.  Unless  the  object  plane  is  large  compared  to  the 
attracted  species  gyroradius,  edge  effects  would  also  be  important.  So,  clearly 
a  much  better,  perhaps  two  or  three  dimensional  treatment,  is  needed,  one  with 
a  more  realistic  sheath  edge  that  integrates  over  all  angles  of  sheath  entry  and 
the  speed  distribution.  Since  such  a  treatment  cannot  be  done  analytically,  it  must 
be  done  numerically.  Such  a  theory  would  resolve  the  specific  questions  raised 


above. 
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